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Abstract 



Close tidal interactions of stars with a central massive black hole (MBH) or with other stars 
\ in the high density cusp around it can affect a significant fraction of the stellar popula- 

tion within the MBH radius of influence. We consider three strong processes that have the 
CO . potential of modifying stellar structure and evolution there. (1) Tidal spin-up by hyperbolic 

star-star encounters. (2) Tidal scattering of stars on the MBH. (3) Tidal heating of inspiraling 
stars — "squeezars" — that were tidally captured by the MBH. We discuss the implications for 
Q-f stellar populations near MBHs and for the growth of MBHs by tidal disruption of stars, and 

Q ■ the possible observational signature of such processes near the MBH in the Galactic Center. 

We compare the event rates of prompt tidal encounters (tidal disruption and tidal scattering) 
and slow inspiral events (squeezars / tidal capture), and find that tidal capture is at least 
an order of magnitude less efficient than prompt disruption. This means that past studies, 
which assigned similar weights to prompt disruption and tidal capture, over-estimated the 
contribution of tidal disruption to the growth of the MBH by at least a factor of two. 

1.1 Introduction 

Strong tidal interactions involving stars are expected to occur frequently near a 
MBH in a galactic center. 

First, the MBH is a mass sink, which drives a flow of stars from the MBH radius of 
influence t% to the center, to replace those it has destroyed. An inevitable consequence of 
this flow is that some stars are deflected into orbits whose periapse r p lies just outside the 
critical radius for destruction. We will focus here on the case where the MBH mass m is 
small enough so that the tidal disruption radius r, ^R+im/M*) 1 ^, where M* and R* are the 
stellar mass and radius, lies outside the event horizon r s , {m< 10 8 M Q for a solar type star). 
Such stars will suffer an extreme tidal impulse, but will not be destroyed, at least not on their 
first peri-passage. There are two possible outcomes: that the star is ultimately disrupted, or 
that it avoids subsequent encounters with the MBH. Both are considered in detail below. 

Second, a variety of formation scenarios predict that MBHs should lie in the center of a 
high density stellar cusp (e.g. Bahcall & Wolf 1977; Young 1980). The diverging stellar den- 
sity implies that there must be some volume around the MBH where close tidal encounters 
occur on timescales significantly shorter than the typical stellar lifetime. Such encounters 
will have a very different outcome from those that occur in globular clusters that do not 
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contain a MBH. In most cases the encounters will not lead to tidal capture. Instead the two 
stars will continue on their separate ways after experiencing a brief strong tidal impulse. 

Extreme tidal interactions, which transfer energy and angular momentum from the orbit 
to the star, can affect its structure and subsequent evolution by heating it, spinning it up, 
mixing it, or ejecting some of its mass. This is interesting in view of the observed presence 
of unusual stellar populations near MBHs: the blue nuclear cluster in the inner ^0.02 pc 
of the GC (Genzel et al. 1997), and around the MBH in M31 (Lauer et al. 1998); evidence 
for anomalously strong rotational dredge-up in an M supergiant near the MBH in the GC 
(Carr, Sellgren & Balachandran 2000), but not in a high density nuclear cluster without a 
MBH (Ramirez et al. 2000); the unusually high concentration of very rare extreme blue He 
supergiants around the Galactic MBH (Krabbe et al. 1991; Najarro et al. 1994). 

The observational consequences of extreme tidal interactions cannot be predicted with 
certainty at this time, although some reasonable conjectures can be made (Alexander & 
livio 2001; Alexander & Morris 2003). Irrespective of this uncertainty, it is clear, as shown 
below, that the amount of tidal energy deposited in the star can reach a significant fraction of 
its binding energy, and that the angular momentum extracted from the orbit can spin a star up 
to a significant fraction of its break-up velocity. It is therefore plausible to assume that the 
effects can be observationally interesting and proceed to explore the dynamical processes 
that give rise to such tidal interactions. Furthermore, tidal disruption and collisional stellar 
mass loss are important channels for supplying mass to a low-mass MBH (Murphy, Cohn & 
Durisen 1991), and so the consequences of extreme tidal processes may provide observable 
links between the properties of the stellar population near the MBH and its evolutionary 
history. 

At a distance of —8 kpc (Reid 1993), the low mass 2.6 x 10 6 M Q MBH in the Galactic 
Center (GC) (Ghez et al. 2000; Schodel et al. 2002) is the nearest and observationally most 
accessible MBH. Although it is heavily reddened (A K ~3 m , Blum et al. 1996), deep high 
resolution astrometric, photometric and spectroscopic IR observations of thousands of stars 
very close to the MBH provide information on their luminosity, effective temperature and 
orbits (e.g. Eckart, Ott & Genzel 1999; Figer et al. 2000; Gezari et al. 2002). Since the 
GC is the obvious first place to look for evidence of extreme tidal interactions, the results 
presented here will be applied to the GC, but it should be emphasized that these physical 
mechanisms are generally relevant for MBH in galactic nuclei. 

1.2 Tidal spin-up by star-star encounters 

Stars in the high Keplerian velocity field in a dense stellar cusp around a MBH will 
suffer numerous hyperbolic tidal encounters over their lifetimes. Although such encounters 
transfer some energy and angular momentum from the orbit to the colliding stars, they rarely 
remove enough energy for tidal capture. This is in marked contrast to the situation in the high 
density cores of globular clusters without a MBH, where the colliding stars are on nearly 
zero-energy orbits and close collisions lead to the formation of tidal binaries. The effects of 
hyperbolic encounters on the stars are mostly transient. The stellar dynamical and thermal 
timescales are very short compared to the mean time between collisions, and so apart from 
some mass-loss in very close collisions, the star is largely unaffected. It is however more 
difficult for the star to shed the excess angular momentum, since magnetic breaking operates 
on timescales of the order of the stellar lifetime (Gray 1992). High rotation is therefore the 
longest lasting dynamical after-effect of a close encounter. Over time, the stellar angular 
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momentum will grow in a random walk fashion due to successive, randomly oriented tidal 
encounters. 

We consider the effect of the tides raised by an impactor star of mass m on a target star 
of mass and radius as the impactor follows an unbound orbit with periapse r p from 
the target star. The tilde symbol will be used below to denote quantities expressed in units 
where G = M* = R*= 1. In these units, 51= 1 is the centrifugal break-up angular frequency, 
Eb = -1 is the stellar binding energy, up to a factor of order unity, and 7, = in 1 / 3 is the MBH 
tidal radius. 

The orbital energy AE and the orbital angular momentum AJ that are transferred from 
the orbit to the star by an impulsive tidal encounter are related by AE = fl p AJ, where tt p is 
the relative angular velocity at periapse. We assume here for simplicity rigid body rotation. 
The change in the stellar angular velocity due to a single parabolic encounter is then given 
to the leading order in the linear multipole expansion by (Press & Teukolsky 1977) 

m p i% 

where / is the stellar moment of inertia (assumed to remain constant), T% is the £ = 2 tidal 
coupling coefficient (calculated numerically for a given stellar structure model), and r] = 
[7^/(1 +m)] 1//2 is the dimensionless transit time of the encounter. The steep ~r~ p dependence 
indicates that most of the contribution comes from very close encounters, where the linear 
expansion no longer holds. The formal divergence at small periapse must be truncated by 
non-linear effects, which have to be investigated numerically. Smoothed Particle Hydrody- 
namics (SPH) simulations (e.g. Fig. I1.2> of grazing and penetrating encounters show that 
as 7 p is decreased, Ail first increases above the value predicted by linear theory, and then 
saturates at the onset of mass loss, as the ejecta carry away the excess angular momentum. 

These results can be incorporated into the linear theory by simple prescriptions, making it 
possible to calculate the mean spin-up of a test star by averaging over all impact parameters, 
orbital energies and impactor masses for a given model of the nuclear stellar cluster (Alexan- 
der & Kumar 1999). Figure ITTTI shows the mean spin-up of a solar mass star after 10 Gyr 
near the MBH in the GC, as function of distance from the MBH, in a high density, n± oc r" 1 5 
cusp of a continuously star forming population that is deduced to exist there (Alexander & 
Sternberg 1999; Alexander 1999). The mean spin-up reaches values as high as ~0.3 within 
0.03 pc, (~60 times higher than is typical in the field) and decreases to ~ 10% within 0.3 pc 
(~0. 1 to 0.2 of ?/,). The spin-up effect falls off only slowly with distance from the black hole 
because the increased tidal coupling in slower collisions at larger distances compensates for 
the decrease in the stellar density. Thus, long-lived main sequence stars with inefficient 
magnetic breaking are expected to rotate at a significant fraction of their centrifugal breakup 
velocity in a large volume of the dense stellar cusp around a MBH . 

1.3 Tidal scattering by the central black hole 

Tidal disruption is an important channel for feeding low-mass MBHs that accrete 
from a low density cusp where collisional mass-loss is low. Numeric models of the growth of 
a MBH in a central cluster suggest that the fraction of the MBH mass that is supplied by tidal 
disruption ranges between /^0.15 (Murphy, Cohn & Durisen 1991) and /^0.65 (Freitag 
& Benz 2002), depending on model assumptions. We will here / = 0.25 as a representative 
value. 
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Fig. 1.1. Tidal spin-up in the Galactic Center (Alexander & Kumar 1999). The 
total mean stellar spin-up of a solar type star after 10 Gyr is plotted as function of 
distance from the MBH. The separate contribution from collisions with other MS 
stars, white dwarfs (WD), neutron stars (NS) and stellar black holes (BH) are also 
shown. 

Dynamical analyses of the 2-body scattering by which stars are deflected into "loss-cone" 
orbits that bring them within 7 t of the MBH (Frank & Rees 1976; Frank 1978; Lightman & 
Shapiro 1977; Magorrian & Tremaine 1999; Syer & Ulmer 1999) show that tidally disrupted 
stars are typically on slightly unbound orbits relative to the MBH, and that they mostly 
originate from the MBH radius of influence, F/, =m/er 2 , where a is the ID velocity dispersion 
far from the MBH. The stellar mass enclosed within ?/, is comparable to in. The cross-section 
for such stars to pass within 7 p of the MBH scales as ~7 r / ,*(Hills 1975; see also 31.51 1. It then 
follows that for every star on an orbit with < r p < 7, that is promptly disrupted, there is a 
star that skirts the tidal disruption zone on an orbit with 7,<7 p < 27, . Such "tidally scattered" 
stars narrowly escape disruption on their first peri-passage after being subjected to extreme 
tidal distortion, spin-up, mixing ad mass-loss (Fig. 11.21 . We will now argue that there is 
also a high probability that these stars will avoid subsequent total disruption, either by being 
deflected off their orbit or by missing the MBH due to its Brownian motion (Alexander & 
Livio 2001). 

First, the scattering timescale is shorter than the dynamical one, and so stars wander in 
and out of the loss-cone several times during one orbital period. After the first peri-passage, 

* This follows from the fact that the stars are on nearly parabolic orbits and that the scattering is isotropic. First, 
since the enclosed stellar mass within roughly equals the MBH mass, a star on a plunging orbit will pass the 
MBH with a velocity slightly above the local escape velocity. Second, is where the loss-cone replenishment 
efficiency peaks sharply, at the transition from the empty loss-cone ("diffusive", or "small angle scattering") 
regime and the full loss-cone ("pinhole", or "large angle scattering") regime. Deflection by large angles relative 
to the loss cone opening angle leads to an isotropic redistribution of the velocity. 
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the stars are on very eccentric orbits with apoastron < 27/,, and so there is a considerable 
chance that they will be scattered again out of the loss-cone before their next close passage, 
and avoid eventual orbital decay and disruption. 

Second, the survival probability is further increased by the Brownian motion of the BH 
relative to the dynamical center of the stellar system. A low mass MBH which evolves in 
an initially constant density core of radius ~r^, is estimated to have Brownian fluctuations 
with an amplitude that is much larger than the tidal radius (Bahcall & Wolf 1976), 

(AT) r h 



55/6 



>1. (1.2) 



The Brownian motion proceeds on the dynamical timescale of the core, which is comparable 
to the orbital period of the tidally scattered stars. The orbits of these stars take them outside 
of 7/,, where they are dynamically affected only by the center of mass of the nucleus, and 
not by the relative shift between the MBH and the stellar mass. Therefore, on re-entry into 
the volume of influence, their orbit will not bring them to the same peri-distance from the 
MBH. 

These order of magnitude arguments are verified by more detailed analysis (see Sll.5> . 
which shows that the survival probability of tidally scattered stars is P s ^0.8-0.9. It then 
follows that the mass fraction of surviving tidally scattered stars within ?/,, which passed 
once within 7 p < 27, of the MBH, is comparable to the mass fraction of the MBH supplied 
by tidal disruption. Depending on the definition of what constitutes an extreme tidal in- 
teraction, there exists a maximal periapse, parameterized by b e = r p /r t , that corresponds to 
sufficiently strong tidal interactions. For example, for a solar type star, b e = 1 .25 corresponds 
to a minimal tidal energy deposition of AE = 0.02 (Eq. II. 4> . Since the angular velocity at 
periapse is f2 p = (2/Z? 3 ) 1 / 2 = 1.01 and the solar moment of inertia 7 = 0.07, this corresponds 
also to a minimal angular spin deposition of A/=0.02 and a minimal spin-up of Af2 = 0.28. 
Overtime, the mass fraction of tidally scattered stars within ^7/, will rise to f(b e —l)P s ~0.05 
(for/ = 0.25, b e = 1.25 and 7^ = 0.8). Tidally scattered stars thus constitute a non-negligible 
fraction of the stellar population in the MBH radius of influence, and will remain there as 
relics of the early stages of the MBH evolution even after its mass grows above the tidal 
disruption limit, possibly detectable by correlations between unusual spectral properties and 
highly eccentric orbits. 

1.4 Squeezars: Tidally powered stars 

A small fraction of the stars that are deflected into orbits with 7 p >7, will be tidally 
caught by the MBH and spiral into an ever tighter orbit as the tides gradually extract orbital 
energy each peri-passage. The orbital energy that a star has to lose to circularize from an 
E = orbit exceeds its own binding energy by orders of magnitude, 

~ m 2 / 3 

£ e = — »1, (1.3) 
2b 

where the periapse is parametrized by b = 7 p /7,. A tidally heated star — a "squeezar" — 
will ultimately be disrupted by expanding beyond its Roche lobe or by radiating above its 
Eddington luminosity. 

The orbital and internal evolution of a tidally heated star in the course of the inspiral 
depends on its initial structure and is coupled to the changes in its mechanical and thermal 
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Fig. 1.2. Snapshots from a Smoothed Particle Hydrodynamics (SPH) simulation 
of a star undergoing an extreme non-disruptive tidal interaction ("tidal scattering") 
as it passes near a massive black hole. Time is measured in units of the stellar 
dynamical timescale. The star passes near the black hole (located outside the frame) 
on a parabolic orbit with a peri-distance 1.5 times larger than the tidal disruption 
distance. Shortly after peri-passage (f = 2) the star appears to be on the verge of 
breaking in two. However, by the end of the simulation, the two fragments coalesce, 
leaving a distorted, mixed and rapidly rotating bound object. 



properties in response to the tidal heating. One approach to the challenging problem of mod- 
eling squeezar evolution is to consider two simplified cases that likely bracket the range of 
possible responses (Alexander & Morris 1993): (1) Surface heating and radiative cooling 
("hot squeezar"), where the tidal oscillations dissipate in a very thin surface layer that ex- 
pands moderately and radiates at a significantly increased effective temperature (McMillan, 
McDermott & Taam 1987). (2) Bulk heating and adiabatic expansion ("cold squeezar"), 
where the tidal oscillations dissipate in the stellar bulk and cause a large, quasi self-similar 
expansion at a constant effective temperature (Podsiadlowski 1996). 

Given these prescriptions, the evolution of the squeezar orbit, size, luminosity and tem- 
perature can be derived from the tidal energy deposition equation 



AE= b / (m»l), (1.4) 
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where R is the expanded stellar radius in terms of the original radius, and the orbital equa- 
tions for the semi-major axis a, the period P and eccentricity <?, 

a = —m/2E, P = 2n^Ja 1 / (1+m) , e=l—7 p /d, (1-5) 

where E is the orbital energy, and Keplerian orbits near the MBH are assumed. 

Figure fOl shows the evolution of a 1M Q hot squeezar, with the tidal coupling coefficient 
calculated for a model of the Sun (Alexander & Kumar 2001). In particular, the squeezar 
evolutionary model relates the inspiral time to to the initial orbital period Po and the periapse 
b. The mean number of squeezars at any given time in the GC is given by n=ioTi, where t o is 
the mean inspiral time, and T,- is the inspiral event rate (estimated in 31. 5> . As shown below, 
n~ 0.1-1 in the GC if the loss-cone is replenished by two body scattering in a spherically 
symmetric system. 

The observational implications of having on average n squeezars near the MBH can be 
expressed by considering the properties of the "leading squeezar" (the one with the shortest 
period). The leading squeezar has, on average, completed f/fo = n/(n+l) of its inspiral. 
Figure fOl also presents the results in terms of the average properties of the leading squeezar, 
as function of n. It is evident that the effects of tidal heating on the leading squeezar can be 
quite pronounced even if n is small, and that the properties of the leading squeezar can be 
much more extreme if n was under-estimated by the neglect of non-spherical effects. 

1.5 Prompt disruption .vs. slow inspiral 

The orbital decay of a tidally heated star is just one case of a dissipative interaction 
that can lead to orbital inspiral. Other possibilities include gravitational wave (GW) emis- 
sion (Hils & Bender 1995; Sigurdsson & Rees 1997; Freitag 2001, 2003) or drag against 
a massive accretion disk (Ostriker 1983; Syer, Clarke & Rees 1991; Vilkoviskij & Czerny 
2002). Unlike prompt disruption or tidal scattering, where the star reaches the MBH directly 
in less than the initial orbital period P , slow inspiral proceeds gradually over a timescale 
to^>Po, which is typically a steeply rising function of the periapse. For the extracted orbital 
energy to power a high luminosity of gravity waves, tidal heat, or mechanical energy in the 
disk, as the case may be, the star has first to decay into a short period orbit. The time avail- 
able for inspiral is limited by two-body collisions similar to those that deflected the star into 
its eccentric orbit in the first place, since they can deflect it again to a wider orbit where the 
dissipation is inefficient. Because to^>Po, this poses a much more severe constraint for an 
inspiraling star than for a promptly disrupted star. 

Novikov et al. (1 992) estimated that tidal capture by a MBH occurs for orbits with b<b c ^ 
3*, which means that stars are scattered into tidal capture orbits and subsequently disrupted 
at a rate that is b c — 1 ~2 times faster than the rate at which they are scattered into prompt tidal 
disruption orbits. These considerations led to the suggestion (Frank & Rees 1976; Novikov 
et al. 1992; Magorrian & Tremaine 1999) that slow tidal inspiral may be at least as important 
as prompt disruption for feeding the MBH and for producing observable tidal flares (Frank 
& Rees 1976). This implies that the already large contribution of prompt tidally disrupted 
stars to the mass budget of a low-mass MBH ( 31.3l l should be further scaled by the ratio of 

* Based on the requirement that the orbital energy extracted by the first peri-passage should decrease the apoapse 
to 2a < ?o, where 7q is the distance from which the star was scattered into the loss-cone. This tidal capture 
criterion does not include timescale considerations. 
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Fig. 1.3. The evolution of a hot squeezar in the GC (Alexander & Morris 2003). 
The 1 Mq star is on an orbit with b = 1.5 and P = 1 .4 x 10 4 yr (f = 4.9 x 10 5 yr). 
The star is disrupted when R = b, at fdismpt = 3.7 x 10 5 yr. At that point the tidal 
luminosity exceeds the intrinsic luminosity by a factor of ~ 640, but the orbit is 
still almost parabolic with eccentricity e= 1 -2.3 x 10~ 4 . The average of properties 
of the leading squeezar as function of the mean number of squeezars can be read 
off the top axis. 

the cross-sections of tidal capture and prompt disruption, b c — 1 . If the relative contribution 
of inspiraling stars were indeed so high, the implications would be far-reaching: stars could 
supply most or even all of the MBH mass, thereby establishing a direct link between in 
and stellar dynamics on a scale of 7/,. However, a small initial periapse does not in itself 
guarantee inspiral and ultimate disruption. The star must also have enough time to complete 
its orbital decay. 

This time constraint can be taken into account correctly by considering only stars that 
are scattered from a volume that is close enough to the MBH so that the inspiral can be 
completed in time (Alexander & Hopman 2003). The inspiral time to increases with T p . 
Therefore, the volume from which a scattered star can inspiral faster than the time it takes 
for two body relaxation to significantly deflect it, decreases with r p . The maximal possible 
periapse corresponds to the point where the available volume shrinks to zero (The truncation 
of the cusp near the MBH by destructive stellar collisions also limits the available volume). 

Table ll . ll lists the inspiral event rate, the mean number of inspiraling stars, and the max- 
imal periapse for inspiral in the GC for three processes: hot and cold 1M Q squeezars and 
for gravitational wave inspiral by O.6M white dwarfs (WD) that comprise 10% of an old 
stellar population. The rates were estimated for a spherical single-mass Keplerian stellar 
cusp, oc f , normalized to contain a total stellar mass 2.6 x 10 6 M Q (= m) within r/, = 1.8 



8 



T. Alexander 



Table 1.1. Inspiral in the Galactic Center. 



Process 


r« (yr 1 ) 


n 


max r p 


Hot squeezar 


3(-6) 


0.2 


2.1 r, 


Cold squeezar 


4(-6) 


0.2 


2.8r, 


WD gravity waves 


2(-7) 


0.04 


25r s 



pc (after Schodel et al. 2002). The predicted prompt disruption rate for this simple model is 
r ; , = 9 x 10~ 5 yr _1 , in general agreement with independent estimates from previous studies, 
T p = 5 x 10" 5 yr" 1 (Syer & Ulmer 1999) and T p ~few x 10" 5 yr" 1 (Alexander 1999). The rate 
of WD inspiral derived here, T, ^2x 10~ 7 yr _1 , is also consistent with the estimates of Sig- 
urdsson & Rees (1997) and Freitag (2003). The survival probability of tidally scattered stars 
can also be calculated by the formalism of Alexander & Hopman (2003), and it is found to 
be close to unity, P s ~ 0.8-0.9, as anticipated by general arguments 

However, we find that the tidal inspiral rate is only ^0.05 of the prompt disruption rate, 
and not the factor 2-3 enhancement due to tidal capture that was assumed by previous stud- 
ies. We conclude that the contribution of tidal capture to the MBH mass budget and to the 
tidal flaring rate from galactic nuclei is negligible compared to prompt disruption. Past stud- 
ies, which assigned similar weights to prompt disruption and tidal capture, over-estimated 
the contribution of tidal disruption to the growth of the MBH by at least a factor of two. 

1.6 Summary 

We have shown that strong tidal interactions of stars with a MBH or with other 
stars in the high density cusp around a MBH can deposit large amounts of orbital energy and 
angular momentum in a significant fraction of the stellar population within a large volume 
around the MBH. We propose that such interactions can alter the evolution and appearance 
of stars in galactic centers and thereby probe the evolution of the MBH and the stellar system 
around it. We explored tidal spin-up by star-star encounters, tidal scattering by a MBH and 
tidal inspiral into a MBH. We showed that tidal capture is inefficient in the presence of two 
body scattering. 
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